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Abstract. The long-time dynamics of the ID contact process suddenly brought 
out of an uncorrelated initial state is studied through a light-cone transfer-matrix 
renormalisation group approach. At criticality, the system undergoes ageing which 
is characterised through the dynamical scaling of the two-times autocorrelation and 
autoresponse functions. The observed non-equality of the ageing exponents a and b 
excludes the possibility of a finite fluctuation-dissipation ratio in the ageing regime. 
The scaling form of the critical autoresponse function is in agreement with the 
prediction of local scale-invariance. 
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1. Introduction 

Consider a statistical system prepared in some initial state. How does it relax 
towards one of its stationary states? This problem was first studied systematically in 
glassy systems where is was observed that the approach towards the thermodynamic 
equilibrium in this kind of system can be very slow (formally the relaxation time 
becomes infinitely large) and depends on the details of the history how the relaxing 
material was treated. While this seemed to preclude any systematic study of such 
systems, it was found empirically P] that the time-dependence of observables can be 
cast into dynamical scaling forms such that universal and reproducible properties of the 
relaxation process emerge. Later, similar effects were also observed to occur in spin and 
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structural glasses (see e.g. j2]) and in non-glassy, e.g. simple ferromagnetic, systems. § 
From a microscopic point of view, the relaxation in these systems proceeds through 
the formation of correlated domains of a time- dependent typical size i{t). This slow 
motion arises because any local discrete spin variable will have at least two distinct 
equilibrium values it could relax into. For a given position r, the locally fluctuating 
variables (e.g. a local magnetic field h{r) in magnetic systems) will select towards which 
of the possible equilibrium states the local order parameter (f){r) will evolve but there 
does remain a competition between the distinct equilibrium states at the macroscopic 
level. For simple ferromagnets and other non-glassy systems, it is generally admitted 
that i{t) ~ t^/^, where z is the dynamical exponent. For systems relaxing towards an 
equilibrium steady-state and which are brought infinitely rapidly from their initial state 
into contact with a thermal bath of temperature T, the value of z depends on whether 
T < Tc or T = Tc, where Tc is the equilibrium critical temperature of the system. It 
is convenient to study these relaxation phenomena through two-time quantities such 
as the two-time auto correlator C{t, s) and the two-time linear autoresponse function 
R{t, s) defined by 

6{<p{t)) 



C(t,s) = (0(t)</.(s)) , R{t, 



(1) 

h=0 



5h{s 

where (f){t) is the time-dependent order parameter and h{s) is the magnetic field 
conjugate to (p. Causality implies that R{t, s) = for t < s. By definition, a system 
is said to undergo ageing, if either C{t, s) or R(t, s) does not merely depend on the 
time difference t = t — s, but on both the observation time t and the waiting time s. 
Ageing occurs for quenches to temperatures either below or at the critical temperature 
but systems in the high-temperature phase with T > Tc do not age. For recent reviews, 
see e.g. El H El El EI • 

Ageing systems may display dynamical scaling in the long-time limit ^ El El IH 
El El EI- Specifically, consider the two-time functions in the ageing regime t 3> tmicro) 
s S> tmicro and r = t — s ^ ^micro, where tmicm is some microscopic time. Then one 
expects the scaling behaviour 

C(t, s) ~ s~'fcit/s) , Rit, s) ~ s-'-'^fnit/s) (2) 

where the scaling functions fc,R{x) have the following asymptotic behaviour for x — *> cxo 

fc{x) ~ x-^^/^ , fnix) ~ x-^«/^ (3) 

Here A^ and \r are called the autocorrelation jSl E] and autoresponse ITUI exponents, 
respectively. The values of the exponents Xc,r and z depend on whether T < Tc or 
T = Tf.. For example, z = 2 for T < T^ and a non-conserved order parameter. In 
general, the exponents Ac and A^j are distinct, but for an infinite-temperature initial 
state it can be shown that Galilei-invariance of the model at zero temperature is a 
sufficient and model-independent criterion for the equality Ac = \r [HI. Indeed, for a 
disordered initial state, Ac = Xr had been taken for granted since a long time, see e.g. 

§ In these studies, it is always assumed that the underlying dynamics satisfies detailed balance. 
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[21 El E]; and for recent reconfirmations in interacting field-theory see ^21 El- In that 
case, one has A^ > d/2 [T3]. The exponents Xc,r are independent of the equihbrium 
exponents and of z (THl . 

For ageing ferromagnetic systems with a non-conserved order parameter, the value 
of the exponent a depends on the properties of the equilibrium system as follows pi)llT7j. 
A system is said to be in class S ii its order-parameter correlator Coq(r) ~ exp(— |r|/^) 
with a finite $, and it is said to be in class L\\ if Ceq{r) ~ ||'|-("'-2+»?)^ where 77 is a 
standard equilibrium critical exponent. Then 

: for class S , ,. 

(4) 
; for class L 

Furthermore, 6 = for T < Tc and b = a ii T = Tc, see e.g. 0111111 • Systems quenched 
to T = Tc are always in class L. 

The distance from equilibrium is conveniently measured through the fluctuation- 
dissipation ratio [IHl dl 

X(t,.):=Ti?(M)(^^) ' (5) 

At equilibrium, the fiuctuation-dissipation theorem states that X(t,s) = 1. Ageing 
systems may also be characterised through the limit fiuctuation-dissipation ratio 

[201 mi El 

Xoo= lim (\imX{t,s)) (6) 




+00 \t— >oo 

Below criticality, one expects Xoo = 0, but if T = Tc, it should be a universal number, 
and this has been confirmed in a large variety of systems in one and two space dimensions 
[211 1221 1231 1121 EH 123 121 • Alternatively, at T = T^ one may fix x = t/s and consider 
X{x) = lims^oo^l^-s, s) and then X^o = ^i^^x^oo X (x) . The order of the limits is 
important, since limf^oo (linis^oo^(^, ■§)) = 1. 

The above discussion of the properties of ageing systems has implicitly assumed that 
the system's dynamics satisfies detailed balance and therefore always relaxes towards 
an equilibrium steady-state even if it may never reach it. Here we wish to investigate 
the relaxation of more general systems where detailed balance is no longer satisfied and 
whose steady-state therefore cannot be in thermodynamic equilibrium. Probably the 
simplest kinetic system fundamentally far from equilibrium is the celebrated contact 
process which has a steady-state transition in the directed percolation universality 
class. ^ 

In section 121 we shall recall the definition of the model and discuss the computation 
of correlators and response functions before we present in section |21 those elements of 
the light-cone transfer-matrix renormalisation group (LCTRMG) which are important 
for our purposes. Calculating two-time observables from the LCTMRG, we study in 
section |3] their time-dependent scaling behaviour. In particular, we shall investigate 

II For example, the kinetic spherical model quenched to T < Tc is in class L. 

% See |27| for an exactly solvable kinetic Ising model with a non-equilibrium steady state. 
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whether there exist scahng forms analogously to eqs. (J2l3p and if so, what the values of 
the ageing exponents a, b, Xc/z, Xr/z are. We shall inquire whether any analogue of a 
finite fluctuation-dissipation ratio might exist. This question is of particular importance 
since in a two-dimensional voter model an effective non-equilibrium temperature can 
indeed be defined through eq. © |2S1- However, the voter model still satisfies detailed 
balance (albeit in a slightly unusual form) and we shall try and see whether the 
procedure proposed in J2S1 niight conceivably be extended to the contact process. A 
more quantitative question concerns the form of the time- dependent scaling functions. 
Indeed, for ageing ferromagnets it has been shown that dynamical scaling as defined 
above can be extended to a richer local scale-invariance [2H]- In particular, the scaling 
forms of the two-time reponse function fniy) 123 ISD] and more recently also of the 
two-time autocorrelation function fc{y) for T < Tc [HI El] have been derived. We 
shall study whether an analogous generalisation is possible in the contact process. A 
complementary paper studies the same model through intensive Monte Carlo simulations 
[S2]- In section El we conclude. 

2. The contact process 

The ID contact process (CP) is defined as follows and might be conceived as a simple 
model describing the propagation of an epidemic disease, see j33l Elj for reviews. 
Consider a chain where each site can be either be empty ('healthy') or else be occupied 
by a particle of a single species A ('infected'). The evolution occurs according to the 
following rules and rates: 

v40,0A^AA ('infection') 
A ^0 ('healing'). 



t O-Oi0-W&.Q Q • O O O O 
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Figure 1. Rules for the eontact proeess on a ID lattiee. Upper line: updating rules 
for a Monte Carlo simulation. Lower line: updating rules for the LCTMRG. 

We shall consider throughout the case of one spatial dimension. Then an empty 
site next to an occupied one is infected with rate A, while infected sites heal at rate 
1, independently of their neighbours. In the upper panel of figure 1 we recall the 
elementary processes for updating a ID lattice with the corresponding rates and this 
readily defines a Monte Carlo algorithm with asynchronous random sequential updates. 
In the lower panel, we present a version of the CP which is more adapted for use with 
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the LCTMRG, see section |21 Besides A, it is also common to parametrise the contact 
process with p := (1 + A)~^ as we shall do in the following: if a lattice site i is occupied, 
the particle is annihilated with probability p, while if the site is empty, a particle is 
created at one of the neighbouring sites of i with probability 1 — p. 
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Figure 2. Density N{t) of the ID contact process in the absorbing phase {p = 0.25), 
at criticality (p = Pc) and in the ordered phase {p = 0.2). Time is measured in units 
of the increment At = 0.02. 



There is a steady-state phase transition between the ordered phase ('infected') at 
small p (large A) and the absorbing phase ('healthy') at large p (small A). To illustrate 
this, consider the average density N{t) := {fiiit)) which is in fact independent of the 
lattice site i and which is computed using the methods described in section |S1 In 
figure 121 N{t) is shown for three values of p. There is a critical value pc = 0.232674(5) 
(A = Ac = 3.29785(8)) [SB] such that N{t) ~ t'^ The exponent 6 ~ 0.16 can be 
read off from the slope in figure in agreement with the literature j^Sl- For p < Pc 
(A > Ac), the density converges exponentially fast towards a positive steady-state value 
A^oo := linit-too -^(^) > 0, while for p > pc {\ < Ac), the density decays exponentially 
fast towards zero. 

Correlations are defined as density-density autocorrelators, both in a disconnected 
and a connected form 



C{t,s) := {ni{t)ni{s)) 

r(t,s) ■.= C{t,s)-N{t)N{s). 

In order to obtain the autoresponse function 

6{ni{t)) 



R{t, s) 



ShAs) 



(9) 



(10) 



/i,=0 



we must introduce an external field hi coupled to the density operator rii at site i. 
The most straightforward possibility is to introduce a spontaneous creation of particles 
— > A on site i with rate hj. 
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Previous experience with R{t, s) comes from magnets, where the system's response 
to an external field is studied via the integrated response since data for R{t, s) are usually 
too noisy. For the integrated response, the most common protocols in magnets are zero- 
field-cooled (ZFC) and thermoremanent magnetisation (TRM), but also intermediate 
protocols have been proposed. However, the analysis of the time-dependent scaling of 
integrated response functions is not always straightforward and the scaling behaviour 
of interest may even be obscured by non-scaling terms or finite-time corrections. This 
occurs in particular for quenches into the ordered phase ^Zj. Another difficulty is 
that the application of a uniform external field h over several time steps introduces 
extra particles into the lattice which brings the steady-state into the active phase. 
For magnets, sophisticated techniques have been proposed [331 EHl EZl I2H] in order 
to avoid such a systematic change of the system's properties and which may involve 
random external fields which change sign in every time step, followed by an average 
over realisations of the randomness or delicate changes in the dynamics. Due to the 
absence of detailed balance in the contact process, these methods are not available here. 

In view of these difficulties, it is a great advantage of the LCTMRG method that 
the autoresponse can be computed directly (see section IHl). This is faster, more accurate 
and conceptually simpler. The limit h ^ can be taken analytically so the problems 
created by a non- vanishing external field are circumvented. 

In calculating C{t,s), r{t,s), and R{t,s), we shall always use a completely filled 
lattice as initial state. This state is completely uncorrelated and can be thought of as 
prepared at p = (A = oo). 

Since the CP does not contain a priori a temperature variable, we might try and 
define an analogy of the fiuctuation-dissipation ratio as follows 

X(t,s):=-4^^^$—. (11) 

^ ' ^ dT{t,s)/ds ^ ' 

3. The LCTMRG method 

Instead of the usual Monte Carlo (MC) simulation of the time-evolution of the contact 
process, see jS21, we use a new variant of the DMRG applied to stochastic transfer 
matrices, the so-called light-cone transfer-matrix density-matrix renormalisation group 
(LCTMRG) [nni- This algorithm, which is an improvement on earlier DMRG approaches 
to calculate the time evolution of stochastic systems [101 E] has the following 
advantages: 

• There is no need for random numbers and ensemble averages since all relevant 
ensembles and correlations (in the sense explained below) are taken into account. 
One LCTMRG run takes a few minutes while 1000 MC runs may take days. The 
resulting correlation functions are very smooth and require no further statistics, 
e.g. in order to compute numerical derivatives with respect to time like dT{t, s)/ds. 

• The transfer matrix enables us to take the thermodynamic limit L ^ oo exactly. 
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However, the LCTMRG is still plagued by numerical instabilities whose exact origin is 
unclear, restricting the calculation to about 1000 time steps. The LCTMRG is not very 
useful for models where each site may have many different states (n ^ 2), or where the 
interaction spans more than two or three sites. 

The dynamics of the one-dimensional stochastic process can be mapped by a 
Trotter-Suzuki checkerboard decomposition onto a two-dimensional classical model: this 
is the geometric interpretation of percolation in two spatial dimensions, directed along 
one of the two axes. The checkerboard is made up of plaquettes r ("local transfer- 
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Figure 3. (a) Trotter-Suzuki decomposition of 2At time steps. The resulting 2D 
lattice consists of local plaqucttc interactions t and is infinitely extended in space 
direction. The dimension of the time direction is finite and the boundary conditions 
are fixed by (l| and |P(0)). (b) Reduction of the 2D lattice to a triangle structure. All 
other plaquettes trivialise, i.e. do not contribute to the state of the top of the triangle. 

After inni- 



matrices" ) encoding the local interaction according to the rules in figure Q 



r 



rir2 



?_2 

{l2r2\e-^'-%n) = rj with /„r,G {0,1} (12) 

h ri 
where h is the local transition-rate matrix from two neighbouring sites hri at time t to 
the same sites at time t + At. The time step At ^ 1 should be chosen sufficiently small. 
We determine the thermodynamic properties of the system by a transfer matrix: 
this ensures that the system is truly infinite in space, while we can follow the short-time 
dynamics for a certain number of time steps. 

Because of probability conservation (Eq. (fT!?jl ) and causality (at each time step, 
only a neighbouring site may be affected by the local interaction), the measurement of 
a local observable riiit) at time step t and site i depends only on the "past light-cone" 
of this site on the classical 2D lattice (|1I]; see figure EJ. 

E/2 Era 

^ = 1. (13) 



V/i,n: > (r):- 
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As the dimension of the exact transfer-matrix grows exponentially with the number 
of time steps, we use the density-matrix renormalisation group (DMRG) idea to decimate 
the state space. The DMRG relies on splitting the system into two strongly correlated 
parts, called the "system" and "environment". Kemper et al. [221 have proposed an 
efficient realisation of the DMRG algorithm applied to corner transfer-matrices. These 
are obtained by diagonal cuts through the checkerboard: the light cone is split into four 
parts diagonally along the future and past light cone of the center point of the triangle 
(see figure E}. For details we refer to [HH] . 




Figure 4. Splitting the light-cone into four corner transfer-matrices (after |39| ) . 



We made a modification to the algorithm in order to be able to compute the 
autocorrelation and autoresponse functions necessary for investigating the ageing 
behaviour. Usually for a local observable rij(t) the expectation value is obtained by 
multiplying the local transfer-matrix at site i and time t with ni{t) before applying 
the initial and final conditions and taking the trace over temporal indices. For the 
two-time correlation function C{t,s), the algorithm has been modified to multiply the 
local transfer matrices r adjacent to site i with rij both at time steps s and t before 
the trace. From C{t,s) the connected autocorrelation T(t,s) is computed via Q, and 
the derivative of the connected autocorrelation function is computed from a symmetric 
difference, i.e. 

dT{t,s) _ T{t,s + At/2) -T{t,s- At/2) 

ds ' At ^ ^ 

which is sufficiently accurate (i.e. independent of At for At = 0.01 . . . 0.05). 

Likewise, when applying an external field hi in order to compute R(t, s), the local 
r adjacent to site i at time step t = s is modified to include particle production at rate 
hi- However, as we are interested in the derivative with respect to the external field, it 
is better to compute this derivative analytically: the Hamiltonian in the presence of an 
external field /i, on site i is 



Hh^ = H + id...,i_i ® hi 



-1 
L 



id 



i+l, 



where H is the stochastic Hamiltonian for the CP, see e.g. |?2l Ho] for reviews. Then 
using the state at time t = s, \P{s)) = e~^'^\P{t = 0)), and the final state (1|, 



Ageing without detailed balance 



R{t,s)= lim (i|e-^(*— ^*') 

hi,At'~*0 



hi At' 



\Pis)) 



li^ ,^^^-Ha-s-At,( iH-H,JAt' + OiiAtr) 
hi,At'^o V hiAt' 



\P{s)) 



(l|e"^(*-^) id...,_i 



-1 
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id 



i+l,. 



in^)) 



where the matrix is written in terms the local basis (0, A) on site i. This has two 
advantages: 

(i) The limit hi —* is taken exactly, thus there is no danger of triggering a phase 
transition by inserting extra particles into the system. 

(ii) No numerical derivative is necessary which would have included the difference of 
two very similar quantities, so this method is numerically more accurate. 
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Figure 5. Comparison of the mean particle-density N{t) found by Monte Carlo (MC) 
and by the LCTMRG for the critical ID contact process, p = Pc- The slope ~S is shown 
the inset. The LCTMRG becomes numerically unstable around t = 1000. Times are 
measured in units of At = 0.02. 



In figure EJ we compare the results for N{t) of the critical contact process in ID 
obtained from Monte Carlo (MC) or the LCTMRG, respectively. First, we observe 
that the LCTMRG data are fairly smooth, as is especially evident when studying the 
exponent S directly. However, we also see that the LCTMRG becomes numerically 
unstable around t = 1000 time steps. This happens because the basis vectors of the 
reduced state space offered during the renormalisation by the LCTMRG method step 
become inadequate: the expectation value of the identity operator (1) is around 1 (as 
it should be) only for the first several hundred time steps, then decreases to below O.L 
However, the onset of instability can in practise always be identified very reliably. The 
reason for this instability is that DMRG works best if system and environment are quite 
strongly entangled, which is not the case here. 
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In this section, we first present the results of our numerical calculations. In the 
subsections 14.41 and 14.51 we shall discuss two important implications of our numerical 
results. 

The MC results in figure E] are obtained using the 3-site update in the first row of 
figure^ with asynchronous dynamics and parallel update, system size L = 10^, 10^ time 
steps, and ensemble averaging over 100 runs. All LCTMRG results are obtained with 
the 2-site update shown in the second row of figure d using the following parameters: 
time step size At = 0.02, number of states in projected state space m = 32, external 
field strength exactly /i — >■ 0, lattice size L = oo (exact thermodynamic limit). 

4.I. Ageing at criticality 

For the critical contact process, we now try and see whether the scaling forms used to 
describe ageing in magnets apply. We consider the scaling forms 



Cit,s) 
T{t,s) 
R{t,s) 



'fcit/s) , 
'fr{t/s) , 
fR{t/s) 



fc{y) 

fr{y) ' 



y 



-\clz 



y 



-Ar/z 



-1-a 



-\r/z 



(15a) 
(15b) 
(15c) 



where the asymptotic behaviour holds for y 

1.0^ 



fR{y) ~ y~ 

cxD. First, we display in figure El the 




Figure 6. Autocorrelation C(t, s) for tlie critical CP for several values of the waiting 
time s as a function of the scaling variable y = t/s. The straight line has the slope 
-0.16 « -(5. 

critical autocorrelator. A nice data collapse is observed with an exponent b = 26 and 
for large values of t/s, we find a power law according to eq. (fT3|) with Xc/z = 6.^ This 
confirms the expectations derived in [22] • 

"*" Here and in what follows, the numbers in brackets give the estimated error in the last given digit (s). 
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Figure 7. Connected autocorrelation T{t, s) for the critical CP, for several values of 
the waiting time s. The slope of the straight line is —1.85. 

Of course, we are interested in the deviations from this as measured by the 
connected autocorrelator r(t, s) which is shown in figured Again, the data collapse 
very well and we obtain the ageing exponents b = 26 and Xr/z = 1.85(10). 
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Figure 8. Autoresponse function i?(i, s) for the critical CP for several values of the 
waiting time s. The full curve labelled LSI is the prediction of local scale-invariance 
fj^{y) = 0.12 • y-i-85(l - l/y)-'^K The inset shows the same data for l<t/s< 3. 

Next, we study the autoresponse function which is displayed in figure |H1 As for the 
correlations, we find a neat dynamical scaling behaviour and for large values oi y = t/s, 
we can read off the exponents 1 + a = 26 and Xr/z = 1.85(10). We shall come back in 
subsection 14 . 51 to a quantitative comparison of the precise functional form of the scaling 
function /i?(?/) with the prediction of local scale-invariance. 

These tests provide evidence in favour of a dynamical scaling of the ageing behaviour 
in the contact process which at least formally is quite analogous to the one found 
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Figure 9. Fluctuation-dissipation ratio (FDR) X{t,s) as defined in eq. Ullfl for the 
critical contact process. For large waiting times s > 100 the curves s~^X{t, s) collapse 
and quickly saturate at w 0.77(1). 



previously in glasses and ferromagnets. However, there are qualitative differences 
between these models and the contact process as becomes explicit when considering 
the exponent relations 



Ar = A 



R 



1 + a 



26 



(16) 



In magnetic systems, the first of these had been simply taken for granted and only 
recently, it has been established (i) that it need not hold for spatially long-ranged 
initial correlations ^2] or in certain glassy systems [H] and (ii) that for phase-ordering 
(hence z = 2) with short-ranged initial correlations local scale-invariance is sufficient to 
actually prove the exponent equality Ar = Xr [TT]. On the other hand, the non-equality 
of the ageing exponents a and 6 at a critical point comes as a surprise. This important 
result can be illustrated in a different way by considering the analogue eq. (fTT|) of the 
fluctuation-dissipation ratio. For the critical CP, we find from figure El for s ^ oo a 
scaling behaviour X{t, s) ~ s- fx{t/s) where the scaling function fx{y) has a finite limit 
value fxioo) ~ 0.77(1). 



4-2. Absorbing phase 

Here, we discuss the behaviour in the subcritical or absorbing phase, with p = 0.6 > Pc 
(A < Ac). The asymptotic behaviour in this phase may be understood by considering 
the extreme case p = 1 first. For p = 1, the particles on different sites are uncorrelated 
and simply decay with a fixed rate. For any fixed site i and with t > s, a moment's 
thought shows that ni(t)ni{s) = riiit) because rij G {0, 1}. Therefore, C(t,s) = N(t). 

This long-time behaviour of the autocorrelation survives in the entire absorbing 
phase, as is illustrated in figures UHl and [TT] for p = 0.6. We see that C(t, s) decays 
exponentially fast but there is no collapse as a function of r = t — s (we point out that 
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Figure 10. Autocorrelation C{t, s) as a function of t = i — s for several values of s in 
the absorbing phase at p = 0.6. The full curve is proportional to exp(— 0.01(t — s)). 
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Figure 11. Autocorrelation C{t, s) as a function of t for the CP with p — 0.6. The 
full curve is proportional to cxp(— O.Olt). 



for intermediate values of s, the data appear to effectively collapse. For p close to Pc, the 
cross-over to the truly asymptotic behaviour may set in very late.). On the other hand, 
when plotted as a function of t, there is a collapse for large values of s, see figure [TT] 
Similarly, r(t, s) = C{t, s) - N{t)N{s) ~ N{t) for s sufficiently large. 

On the other hand, for the response function we find in figure ^1 a rapid collapse 
in terms of the time difference t = t — s for larger waiting times s > 100 and time- 
translation invariance is recovered, as it should be. 



4-3. Active phase 

Finally, we discuss the long-time behaviour of the two-time observables in the active 
(percolating/ordered) phase, where p = OA < pc {\ = 9 > Xc). In distinction with 
magnetic systems, where ageing also occurs for quenches into the ordered phase, in the 
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Figure 12. Autoresponse R{t, s) as a function of t = i — s for several values of s in 
the absorbing phase at p = 0.6. The full curve is proportional to exp(— 0.01(t — sj). 



contact process time-translation invariance is rapidly recovered and no ageing occurs. 
We illustrate this in figure ^1 for the connected autocorrelator and in figure ^^ for the 
response function. In both cases, when plotted against t = t — s, we observe a collapse 
for waiting times s > 50 and time-translation invariance is recovered. 
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Figure 13. Connected autocorrelation T{t, s) in the active phase of the contact process 
(with p = 0.1). The straight line is proportional to exp(~0.05(i — s)). 



4-4 ■ Effective temperature 

One of the main peculiarities of equilibrium states with respect to more general steady- 
states is that their probability distribution has the form of a Boltzmann weight which 
may be characterised in terms of a temperature. It therefore is a natural question 
whether for more general steady-states a non- equilibrium temperature might be defined. 
Indeed, such an attempt has been presented recently by Sastre et al. [221 • They 
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Figure 14. Autoresponse R{t, s) in the active phase of the contact process (with 
p = 0.1). The straight line is proportional to exp(— 0.05(t — s)). 



considered the fluctuation-dissipation ratio eq. (0) and observed that in the hmit s ^ oo 
and t — s ^ 0, one should have an equilibrium-like regime, hence X{t,s) -^ 1. From 
this observation, they define a dynamical temperature by 
1 ,. / .. R{t,s) 



-dyn 



lim 

t— ►oo 



lim 

t-s^odC{t,s)/ds 



(17) 



By explicit calculation, they confirm that in the 2D critical voter model this limit exists, 
has a non-trivial value and is universal j^Sj- However, the basic assumption of the idea 
of Sastre et al. has been critically reexamined by Mayer and Sollich J3S! who construct 
in the ID Glauber-Ising model undergoing coarsening a defect-pair observable such that 
the fiuctuation-dissipation ratio X{t, s) 7^ 1 in the short time- regime (in particular they 
show \im.s^ooX{s, s) = 3/4). 

We now ask whether the definition ()17|) J2S] could be extended to the contact 
process. Obviously, for the contact process the connected correlator r(t, s) must be 
used in eq. (fTTjl . Before we do so, it may be useful to recall some constraints on the 
ageing behaviour which hold for relaxation towards equilibrium. Equilibrium states are 
steady-states (hence time-translation invariant) which satisfy in addition the fiuctuation- 
dissipation theorem (thus X(t, s) = 1). This last condition may also be replaced by the 
Onsager symmetry condition {A{t)B{s))c(i ~ (-^(^)^('^))cq for ^"^J two observables A 
and B; which together with time-translation invariance is enough to reproduce the 
fiuctuation-dissipation theorem P^IHj. 

We are interested in equilibrium critical dynamics. Combining the scaling forms of 
section^ and time-translation invariance, we expect 

C(t,s)~(t-s)~' , /?(t,s)~(t-s)~'"" (18) 

Combining this with the FDT, it follows a = b. Therefore, the equality a = b is a. 
necessary condition that the quasi- stationary state found when t — s <^ 1 is indeed a 
quasi-equilibrium one. 
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We can apply this argument to the quasi-steady state of the contact process 
characterised by the connected autocorrelator T(t,s) (instead of C{t,s)) and the 
autoresponse R(t,s). In subsection 14. ![ we have seen that dynamical scaling eq. ()15j] 
holds with the ageing exponents satisfying 1 + a = b = 26. Consequently, since a 
necessary condition for the existence of a quasi-equilibrium regime is not satisfied in the 
contact process, the definition (fTTj) of a non-equilibrium temperature cannot be extended 
to this universality class. It is remarkable to see that the presence or absence of detailed 
balance leads to very different response behaviour, once the dynamical scaling-regime is 
reached. We conclude that the definition (fT7|) proposed in [23] is likely to reflect peculiar 
properties of the critical voter model rather than being generic. 
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Figure 15. Ratio R{t, s)/T{t, s) for the critical CP. For s > 25 the curves appear to 
collapse to a scaling function which in turn saturates at a finite value for t/s > 10. 



Indeed, our data suggest that not the 'fluctuation-dissipation ratio' X{t, s) as 
defined in eq. flllj) . but rather the ratio 

R{t,s) 



E{t,s) 



T{t,s) 



(19) 



should become a scale-invariant function of t/s in the ageing regime at criticality, viz. 
S = S(t/s). This appears to be the case as is shown in figure IT31 and for t/s -^ oo 
we read off the limit value Sqo — 1.15(5). In analogy of the universal limit fluctuation- 
dissipation ratio Xoo |2B i^ systems with detailed balance, it would be interesting to see 
whether or not Sqo is universal (for example by going beyond the mean-field treatment 
|32] of Reggeon field-theory) . 



4.5. Local scale-invariance 

Having seen that at criticality, the ageing of the contact process satisfies dynamical 
scaling, it is natural to inquire whether the recently proposed extension [2H] of dynamical 
scaling to a space-time-dependent, i.e. local kind of scale-invariance in magnets might 
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also apply to the model at hand. The central assumption of that theory is that the linear 
response function formed from so-called quasi-primary fields transforms covariantly 
under the action of the group of local scale-transformations. We propose to try and see 
whether the form of the linear response in the critical contact process may be understood 
this way. 

Indeed, for any given value of z, infinitesimal local scale transformations with a 
space-time-dependent rescaling factor 1 + e{t, r) can be constructed. In particular, the 
following explicit expression for the response function is obtained from the condition that 
R{t, s) transforms covariantly under the action of local scale transformations 



l+a-Xn/z 



R{t,s) = roi-j (t-s)-'-'^ (20) 

and where tq is a normalisation constant. This prediction has been confirmed in several 
models with a dynamics given by a master equation, notably the kinetic Ising model 
with Glauber dynamics, both in the bulk [23120] as well as close to a free surface [IB], the 
kinetic XY model with a non-conserved order parameter |Tri [TT| and for the Hilhorst-van 
Leeuwen model ^^1- Local scale-invariance has also been confirmed for several variants 
of the exactly solvable spherical model and the free random walk [221 1^ ESl lllj. 

All these tests are based on a Master equation or on a linear Langevin equation 
which reduces to a free field-theory. On the other hand, field-theoretical calculations 
of ageing ferromagnetic systems based on non-linear Langevin equations find small 
corrections to eq. (|2n|) |T2lEj. However, it is not completely clear whether a description 
of ageing in terms of a Master equation and in terms of a Langevin equation are 
completely equivalent. A counterexample is provided by the ID kinetic Ising model 
with Glauber dynamics at zero temperature where the universal exponent Ac = 1 is 
known exactly fIU\ 1^ . On the other hand, the Langevin equation usually believed 
to be equivalent to this model, namely the time-dependent Ginzburg-Landau equation, 
leads to the exact result Ac = 0.6006 . . . |^ . 

Comparing eq. (PUI) with the LCTMRG data for the ID critical contact process, 
we find from figure |H1 a perfect agreement almost down to t/s = 1, where values for the 
exponents a and Xr/z determined previously were used.* The contact process hence 
provides the first example of a model satisfying local scale-invariance which has a non- 
equilibrium steady-state. 

5. Conclusions 

We have studied the question if an analogue of the ageing phenomenon well-known in 
magnetic systems also occurs in systems which relax to a non-equilibrium steady-state. 

* Recently, it has been established that in those ageing systems which undergo cluster dilution rather 
than domain growth, there is an universal early-time regime where dynamical scaling does not hold 
|51| . On the other hand, from Monte Carlo simulations we know that in the contact process ageing 
proceeds via cluster dissolution [321 Fig. 2]. Therefore a small deviation of the numerical data from 
the t/s scaling for t/s w 1 is to be expected. 
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d 


a 


b 


\c/z 


Xr/z 


\r/z 


method 


1 


-0.68(5) 
-0.57(10) 


0.32(5) 
0.319 


0.16(1) 
0.159 


1.85(10) 
1.9(1) 


1.85(10) 
1.9(1) 


LCTMRG 
Monte Carlo 


2 


0.3(1) 


0.901(2) 


0.450 


2.8(3) 


2.75(10) 


Monte Carlo 


>4 


i-1 


2 


1 


- 


f + 2 


mean-field 



18 



Table 1. Nonequilibrium exponents as defined in (|15|l for the ageing of the critical 
contact process in several space dimensions d, according to Monte Carlo |321 or 
LCTMRG (this work) calculations. The exponents following from mean-field theory 
[^21 are also included. 



Our case study of the one-dimensional contact process has allowed us to answer this 
general question affirmatively, at least in situations when the steady-state is critical and 
might be hence viewed as being formed from two coalescing steady-states. On the other 
hand, in the active phase, where only a single stable steady-state exists, time-translation 
invariance is rapidly recovered and no ageing occurs. At criticality, the ageing behaviour 
can be described in terms of dynamical scaling, see eq. (fTH|l . and we collect the values 
of the exponents in tabled It is satisfying to see that in ID we find a good agreement 
with the results of a Monte Carlo study jH21 on the same model, jj While we find evidence 
for the exponent equality Ar = Xr in close analogy to magnetic systems quenched to 
criticality from a fully disordered state, the ageing exponents a and b are different and 
we conjecture 

l + a = b = 26 (21) 

where (5 is a well-known non-equilibrium exponent. In particular, this implies that 
the quasistationary regime of the contact process is already out of equilibrium and 
the recent attempt j23] to define a genuine non-equilibrium temperature does not go 
through. Instead of the fluctuation-dissipation theorem valid for equilibrium systems, 
we have found evidence that in the contact process rather the ratio R{t, s)/r{t, s) should 
converge to a finite value in the limit of widely separated times (see figure ITH|l . 

On the other hand, we have seen that the scaling form of the response function is in 
agreement with local scale-invar iance. We point out that this confirmation is obtained 
in a formulation based on the master equation and not in a field-theoretical setting 
based on a Langevin equation. All existing confirmations of local scale-invariance, see 
[Z|, have either been obtained in this setup or else come from models which reduce to 
free field-theories. 
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tt This agreement is not completely trivial, since the initial densities in these two studies are different. 
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